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The Hilbert spaces of particle states are constructed in such a way that the neutrinos are combined 
in a multiplet with its components being considered as different quantum states of a single particle. 
The same is done for the charged leptons and the down- and up-type quarks. In a theory based 
on the Lagrangian of the fermion sector of the Standard Model modified in accordance with this 
construction, the phenomenon of neutrino oscillations arises as a direct consequence of the general 
principles of quantum field theory. By the example of the pion decay it is shown that the states 
of the neutrino produced in the decay process can be described by a superposition of states with 
different masses and identical canonical momenta with very high accuracy. 


INTRODUCTION 


The Standard Model of electroweak interactions based on the non-Abelian gauge symmetry of the interactions 
[I], the generation of particle masses due to the spontaneous symmetry breaking mechanism 0 , and the ideology of 
particle generation mixing 1,1 is universally recognized. Its predictions obtained in the framework of perturbation 
theory are in a very good agreement with the experimental data. There is no serious reason, at least at the energies 
available at present, for its main propositions to be revised. 

However, in describing such an important and firmly experimentally established phenomenon as neutrino oscillations 
(see |5|), an essentially phenomenological theory based on the pioneer works by B. Pontecorvo and Z. Maki et al. 
j m| is used. Its statements are given in various review articles and monographs (see, e.g., 0, i, 0, 0, 0, 
0,111, 0). The primary assumption of this theory is that the neutrinos are massive and there are three neutrino 
types with different masses. It is also postulated that the neutrinos produced in reactions are in the states which are 
superpositions of states with fixed masses, the so-called mass states. These states form the so-called flavor basis. The 
transformation to this basis from the mass basis is given by a unitary mixing matrix. Initially, the mass basis elements 
are described by plane waves with the same (three-dimensional) momentum. The time evolution of the flavor states 
is described by a solution of the corresponding Cauchy problem. 

The conclusions of this theory hold true when the neutrino energy is large compared to their masses <C 1. 

Neutrino masses are extremely small. Due to this fact, neutrino energies, which satisfy this condition are not too high. 
Therefore, there is an obvious logical contradiction. On the one hand, neutrino oscillations are observed at energies 
that are typical for the Standard Model. On the other hand, oscillations suggest the possibility of transformations of 
free fermions with equal electroweak quantum numbers to each other, contrary to the the Standard Model. 

Numerous attempts based on the quantum Held theory have been made to justify the conclusions of the phenomeno¬ 
logical model. It is worth mentioning the description of the states of neutrinos in the form of wave packets (see, e.g., 
14]). Another approach involves a description of the birth of the neutrino and its detection as a single process (see, 
e.g., [0, [0). The review of the models based on the quantum theory of fields is given in the paper [18| . The results 
obtained partially explain the oscillation mechanism. However, the theory in the mathematical sense is not closed and 
has a number of paradoxical properties that are clearly formulated in [19j . 

All the contradictions of the theory could be solved, if we could combine the neutrinos in multiplet with its 
components being considered as different quantum states of a single particle. This could be done for the charged 
leptons and the down- and up-type quarks as well. With this purpose it is necessary to construct the Fock spaces for 
flavor states of the particles. 

Attempts to construct the Fock space for flavor states of neutrinos have been repeated many times [0, [21], [0 . 
However, either commutation relations for operators of creation and annihilation differ from the canonical type, or 
the vacuum of the theory was clearly dependent on the time. These investigations resulted in the conclusion that it is 
impossible to construct the Fock space for the flavor states using the conventional approach. In other words, the mass 
and flavor states cannot be connected by a unitary transformation due to introduction of the mixing matrix only. 
Note that, such a possibility was taken for granted in the phenomenological theory. The discussion of the problem is 
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contained in the summary report |23[ . 

Let us consider this issue in detail. In relativistic quantum field theory a particle is usually associated with an 
irreducible representation of the Poincare group [H| , (25j . The eigenvalue of the Casimir operator constructed from 
the translation operators squared is identified with the observed mass of the particle. The existence of states that 
are superpositions of states with different masses contradicts the relativistic invariance of the theory, if the canonical 
momentum operator is identified with the translation operator. This obviously follows from the fact that the metric 
is defined on the hyperboloid in the momentum space determined by the value of the particle mass. 

To overcome this difficulty it seems natural to associate some set of particles (a multiplet) with an irreducible 
representation of a wider group. A number of theorems f26[ indicates that the only reasonable extension of the 
symmetry group of the theory is the direct product of the Poincare group and a group of internal symmetry. This 
immediately limits the set of the models under consideration. However, if we assume that the masses of particles in a 
multiplet can be set by hands, as it is done in the Standard Model, this circumstance is not essentially important. It 
is essential that in this point we encounter the Jost theorem 


27], |28[. It seems that due to this theorem the result of 


extending the symmetry group of the theory will be trivial. Since the irreducible unitary representations of the direct 
product of groups are unitary equivalent to the tensor product of the representations of these groups with the same 
values of the invariants [29[ . the masses of the components of the multiplets will be equal. Just this circumstance is 
the most important reason of the failure of constructing the Fock space for flavor states. 

In our opinion, it is possible to circumvent this obstacle. The basic idea of the proposed approach is as follows. 
The conclusion mentioned above is not correct if there is more than one multiplet in the theory, and these multiplets 
interact with each other. As a matter of fact, in the general case the representation spaces for two multiplets cannot be 
simultaneously reduced to the direct sum of representation spaces of the Poincare group. Even if the representation 
of the internal symmetry group is considered finite that seems quite natural for a finite set of particle masses in 
the multiplet, the representations for both multiplets would be unitary equivalent, but not identical. The values of 
the invariants related to the translation operators squared will be the same for all component of the multiplet (for 
different multiplets they may vary). However, there is no need to treat them as the observed squared masses of the 
particles. The masses will be the eigenvalues of an operator in the enveloping algebra of the Lie algebra of the internal 
symmetry group, which naturally commutes with the elements of the Lie algebra of the Poincare group. Therefore, 
these representations, being unitary equivalent are not “physically equivalent”. Such consideration is conformed to the 
spirit of the Standard Model, where the masses of all particles are generated due to the phenomenon of spontaneous 
symmetry breaking and are proportional to the vacuum expectation value of the Higgs field. 

The purpose of this work is to realize these ideas for developing a modification of the Standard Model, which is 
suitable to describe the oscillations of the particles. To do this, in Section U we construct wave functions of multiplets 
of fermions with spin 1/2, which form a spaces of irreducible representations of the direct product of the Poincare 
group and the group SU( 3). First, we consider the most simple version of such a construction that is the tensor 
product of the Dirac representation of the Poincare group and the fundamental representation of 517(3), which we 
designate as ■ In this space the wave functions describe multiplets with components possessing equal masses 

m. 

Then, on the example of neutrino multiplet, we find unitary intertwining operator mapping 'H m ^\/2 to a new space 
In this space the wave functions describe multiplets with components possessing different masses that are 
the eigenvalues of the operator, which is defined by a mass matrix of the multiplet. These components describe the 
so called mass states of the field. The mass states generate a basis in but the choice of this basis is not 

unique. We find all the plane-wave bases, which can describe the mass states. For any choice of the mass basis, linear 
combinations of its elements are bases in Due to unitary equivalence of representations in 'H rn ,i /2 and in 


^■rn 1 / 2 ' eigenvalues of the Casimir operator constructed from the translation operators squared on the elements 

of 'H^\j 2 are equal to m 2 . Therefore, the parameter m is not the observed mass, it sets the scale of the multiplet 
masses. As a matter of fact, the parameter m has the meaning of a vacuum expectation value of the Higgs field. 

In Section|TTJ using the results of Section|T| we introduce a Lagrangian generalizing the fermion sector of the Standard 
Model in such a way that field functions of the individual particles are replaced by field functions of the multiplets. 
In the model under consideration the action defined by the Lagrangian of free fields is invariant with respect to the 
direct product of the Poincare group and 517(3). Using Noether’s theorem we find the integrals of motion for the 
fields and carry out the quantization procedure. We make sure that creation and annihilation operators satisfy the 
canonical commutation relations. However, these operators possess an additional discrete quantum number that is 
associated with the mass of the state. The results of the quantization make it possible to treat both the mass states 
of the multiplet and their superpositions as the quantum states of a single particle and to construct the appropriate 
Fock space. Therefore, the transition probabilities can be calculated in the framework of perturbation theory. 

In a theory based on the Lagrangian of the fermion sector of the Standard Model modified in accordance with this 
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construction, the phenomenon of particle oscillations arises. In the Standard Model the matrices of Yukawa couplings 
are introduced to describe the interaction of fermions with the Higgs field. In the framework of Kobayashi -Maskawa 
formalism [5j], a diagonalization of these matrices is necessary to generate masses of fermions. This transformation 
leads to a mixing matrix in the charge-current interaction Lagrangian. This procedure presupposes that all fermions 
describing by the model (i.e. the neutrinos, the charged leptons and the down- and up-type quarks) are fundamental 
particles. Therefore, although the mixing of generations occurs, the direct transitions between particles with the same 
electroweak quantum numbers, i.e. oscillations, are impossible. 

In the proposed model mixing occurs at the level of perturbation theory. So there is no need to carry out the 
diagonalization procedure and to introduce the mixing matrix into the interaction Lagrangian explicitly. It is obvious 
that the transition probabilities, obtained in this way, will oscillate. However, at least the transition probabilities 
obtained in a tree approximation will be completely coincide with the predictions of the conventional approach at 
distances from the source that are small compared with the lengths of oscillations. 

In the following sections we apply the results to the description of neutrino oscillations. In Section IIIII we check 
that the formulas describing the neutrino oscillations are independent of the type of the mass states, i.e. of the choice 
of the plane-wave basis in In Section m we calculate the probability of the pion decay and check that the 

probability of creation of one or another flavor state essentially depends on the types of the mass states of which the 
flavor states are constructed. A single flavor state is produced only when it is composed of the mass states with the 
same canonical momentum. Thus, we come to the conclusion that the phenomenological theory of oscillations based 
on the ideas of Pontecorvo is a very good approximation for quantum field theory. 


I. SPACES OF WAVE FUNCTIONS 

Let us assume that the neutrinos, the charged leptons, the down- and up-type quarks are elements of different 
multiplets. For each multiplet we consider a space 'H rrii 1 / 2 , which is the direct sum of the Dirac representation spaces 
for the Poincare group with fixed (positive) frequency and equal values of the Casimir operator formed from the 
canonical momentum operators squared. Since the experimental data indicate that there are three generations of 
particles, 

Km, 1/2 = «£ 1/2 ® n™ 1/2 © i /2 , (1.1) 

where the subspaces ’H'di 2 correspond to the components of the multiplet. 

In the coordinate representation a basis of TLm, 1/2 can be defined as 

%,cA x ) = i>P,d x ) e(i) • (1-2) 

Here ip p ^{x) are plane waves 

%,d x ) = yp u ^ e-i(px) ’ 

where p° = \/p 2 + m 2 , and the spinors u P x satisfy the equation 

(7 *Pn ~ m)u p ,c = 0. (1.4) 

This spinors are normalized by the condition 

%,C u p,C'= 2m ^C,C> UpX = u lx"f° ■ (1-5) 

The indices £ = ±1 determine a particle polarization. 

We can consider the space of the representation with negative frequency, too. A basis in this space can be introduced 
as 

P ,cA x ) = i>-p,d x ) e(0 > i>-p,d x ) = y= M - P ,c el(px) i ( 1 6) 

+ m ) u -pX = W-p,c u -p,C' = -2?nS^x'- 

The explicit formulas for the spinors u p ^, u p .(-, and the spinors U- P ^,U- P ^, see, e.g., in [3lj . In what follows we will 
work in the spaces of wave functions with the positive frequency. The conversions for the spaces of wave functions 
with the negative frequency are similar. 





Consider the vectors which, for definiteness, can be chosen in the form 


e 


(!) = 





(1.7) 


as a basis of a three-dimensional vector space over the field of complex numbers. Introduce a scalar product in "H mi 1/2 

with summation over the coordinates of the vectors e^. Then one can consider 'H rrit i /2 as the space of an irreducible 
unitary representation of the direct product of the Poincare group and SU( 3). This space is constructed as the direct 
product of the Dirac representation of the Poincare group and the fundamental representation of the group SU (3). 

The explicit form of the Lie algebra elements in the present case is obvious. We have the standard realization for 
the generators of the Poincare group 


dx.^ (x, t) #i(x,i), 


( 1 . 8 ) 


Pp = i<9 M I, M )iV = + (1.9) 

where I is the 3x3 identity matrix. The Hermitian generators Xk of the fundamental representation of the SU (3) 
group are defined by the Gell-Mann matrices 

X 1 = ( e W <g> e < 2 )) + (>)<g>e«), 

X-2 = i (e^ <E> — i ® e^ 2 ^ , 

X 3 = ® ^e ^ ® , 

X 4 = fe (1) ®e (3) )+ fe (3) ® e (1) ) , 

V / V / Cl 10) 

X 5 = i (e< 3 > ® e«) - i (e« ® e< 3 >) , 

X 6 = ^e^ ® e^ 3 ^ + ^ 3 ) ® e (2 ^ , 

X’j = i ® e^ 2 ^ — i ® , 

X 8 = (e«®e«)+ ( e ( 2 ) ® - 2 (e® ® e< 3 >) . 

The irreducibility condition for this representation is the matrix Dirac equation (see, e.g., M) 

(iy^c^I — ml)&(x) = 0. (1.11) 

Obviously, in this representation all the components of the multiplet have equal masses. 

Let us turn to another representations of the extended symmetry group of the theory for each multiplet. The idea 
of the proposed transformation is based on the fact that the derivation algebra of the Poincare algebra contains not 
only the operators P 11 , M llu , but also the generator of dilatation D (see, e.g., ,3t|) 

[P M , D\ = P M , [AP", L>] = 0. (1.12) 

Because of this, it is possible to construct an external automorphism of the direct product of the Poincare group and 
SU( 3). For the representation in the space the automorphism leads to scaling transformations of coordinates. 

The transformations can vary for different subspaces. As a consequence, it will allow to consider the multiplets with 
different masses of the components. 

To be specific, we assume that we work with the neutrino multiplet. Introduce new basis vectors n® with the help 
of a unitary matrix acting on 
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The new representation space / H < ^\/ 2 will a l so be defined as the direct sum of the spaces of the Dirac representations 
of the Poincare group 


-T/C^) _ 0 /(^ 1 ) /T\ ^ 

n m,l/2 ~ n m, 1/2 ® 1/2 ® n m, 1/2' 


We define a basis in T~L 


W 

m, 1/2 


in the form 


where ip^ 1 '} (x) are the plane waves derived from (11.31) by a dilatation of the coordinates 


= J ^ Uq ^ e 

Here q° = 7/q 2 + m 2 , and //; are positive real numbers. The spinors u q ^ satisfy the equation 

( 7 - m ) u qX = °> 

and are normalized by the condition 

/(") 


Mi 


3/2 




(1.14) 


(1.15) 


(1.16) 


(1.17) 


(1.18) 


Evidently, (the elements of the space x / 2 ) anc ^ ^( a; ) (tbe elements of the space ' H m> i/ 2 ) are connected 


by the unitary (with respect to the scalar product (11.81) ') transformation 

^ v \x) = JC^^(x) = J K(x, y) ^{y) dy. 
Its kernel is determined by the formula 


K ( X ’V) = (^jlE E j ® KxM- 


The explicit form of the dilatation operator is 


D = x a d a + 3/2. 


(1.19) 


( 1 . 20 ) 


( 1 . 21 ) 


Therefore, the intertwining operator for the representations in the spaces H^ 1 / 2 and realizing the transfor¬ 

mation CLU is 


where 


KS v) =Y j V {1) 

1= 1 


(Dp) = exp (In [ii(x a d a + 3/2)). 


( 1 . 22 ) 


(1.23) 


A direct calculation yields the explicit form of the elements of the Lie algebra of the considered representation. 
The action of the generators X ( k v> = /C^XfcXA'b -1 of the SU( 3) group on the basis (11.151) is reduced to obvious 
permutations of its elements 


vMtp-W i x \ _ 1 p( v ) ( x \ (x) = ^ v } (x) (x) = 0 

• A 1 W qX^X X > Q, W q. W q, C,tW X '> ^1 W C,A*3 W 


(1.24) 


and so on. The diagonal operators are determined by numerical matrices. The explicit form of these 

operators see in appendix l~A~l 

The generators of the Lorentz group do not change the form 

= /C^M^/CH - 1 = \{(x^d v - x„d M ) + ( 7 ^ 7 ,, - 7 „ 7 m )/ 4 ^I, 


(1.25) 
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since the dilatation operator commutes with the generators of rotations and boosts (see (11.121) '). This fact is quite 
natural, because the results of observations, in particular, of the study of particle oscillations cannot depend on the 
choice of the inertial reference frame used for measurement. 

For the translation generators, we have 




where 





(1.26) 


(1.27) 


That is /q 1 are the eigenvalues of the matrix I'b"), and n^ are its eigenvectors normalized by the conditions (the 
asterisk denotes the complex conjugation) 

X>?>sf=4„ = < L28 > 

i=1 1=1 

The operators 

= n m <g) h (l) (1.29) 

are orthogonal projectors 

( L3 °) 

2=1,2,3 


The Dirac equation, i.e. the irreducibility condition for this representation, is now written as 

(ry^lN^ - ml) &M(x) = 0. 

Since the parameters /q are chosen nonzero then IN*-") is non-degenerate and the inverse matrix exists 

1=1 

If we multiply ()1 .31 j) by (|1 .32[) then the Dirac equation takes a more familiar form 

(i7 M <9 M I-MMW'')(, r ) =o, 


(1.31) 


(1.32) 


(1.33) 


where 


3 

M *- 1 ^ = to ( 1 N ( ' ! ')) _1 = 1711 = ^ im " ( 1 - 34 ) 

i=i 

The sets of solutions (11.311) and (11.331) coincide. Further, we will refer only to Eq. (11.331) . 

We emphasize once again that the eigenvalues of the Casimir operators constructed from the generators (11.251) . 
(I1.26D , and (IA1I) take the same values on the solutions of these equations and on elements of the initial representation 
space 'Hm, i/ 2 - In particular, = m 2 <F^^(a;), but now the parameter m is not the observed mass, it 

rather sets the scale of the multiplet masses. In other words, it represents the bare mass of the multiplet. The observed 
masses are determined by the action of the canonical momentum operator squared on the functions ( x) and 

are equal to mi = Him for i-component of the multiplet. 

Let us discuss this point in detail. The general solution of (11.331) can be expanded in the basis functions that are the 
eigenfunctions of a complete set of operators commuting with the operator of the equation. Naturally, the choice of the 
complete set is ambiguous. For the standard Dirac equation, if we restrict ourselves to the case of plane-wave solutions, 
the complete set includes the space components of the momentum (i. e. the generators of the space translations) and 
the operator of the spin projection. This operator is constructed in the common way from the components of the 
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Pauli Lubanski-Bargmann vector. The ambiguity of the complete set composition in that case depends on the choice 
of the spin operator 32]. In contrast to the standard Dirac equation, a complete set for Eq. (11.331) contains five 
operators. This circumstance makes the situation more complicated. 

We now find the complete sets of operators for which the basis of the solution space of Eq. (11.331) consists of plane 
waves. Let P*-") be a non-degenerate Hermitian matrix which commutes with This matrix can be written as 


pM = 


^ A,P^ 


i=i 


(0 


A i £ 0. 


(1.35) 


If we assume that the parameters A/ are pairwise distinct, then, in the space of the solutions of Eq. (11.331) . this matrix 
defines the action of the operator which is a linear combination of (the third projection of the T-spin), 

(the hypercharge), and the Casimir operators of SU(3). (Since the representation is irreducible, then the action of the 
Casimir operators is given by the matrices that are multiple of the identity matrix). Obviously, M ^ can be chosen 
as one of the operators of the complete set. This operator is used to isolate the orthogonal subspaces determined by 
the projectors Pj]/ 1 . Therefore, the numerical values of the parameters A; can be arbitrary. It is natural to set A i = mi, 

that is, to assume that the action of the operator M 1 -"' 1 is defined by the matrix M("). 

For a basis to be a plane-wave one, the complete set must include three operators with continuous spectrum. If 
we choose the space translation generators, that is, the spatial components of the 4-vector P^\ as such operators, 
then any standard spin projector multiplied by IN^ will determine the spin projection. In this case the complete 
orthonormal system of solutions of Eq. (11.331) corresponding to a fixed (positive) frequency coincides with the basis 


& 


m( x ) (see Eq. (11. 15(1 b which was considered previously. 
However, we can take the spatial components of the 4-vector 


(1.36) 

to play the role of the operators with continuous spectrum, where p( y ) is any matrix which satisfies (11.351) . In particular, 
we can set lb") = (USb^ ) _1 . Then the basis functions are the eigenfunctions of the spatial components of the canonical 
momentum operator id and 


^ v l (x ) = ib^l (x)n 
where mj (x) are the plane waves describing particles with masses mi 




(x) 




-i(p*) 

V c e 


(1.37) 


(1.38) 


(1.39) 


(1.40) 


Here p° = -^/p 2 + mf, the spinors u^' 1 satisfy the equation 

( 7 - m)u { ™ c l) = 0 , 

and are normalized by the condition 

= 2m 

Thus, the matrix can be interpreted as the mass matrix of the neutrino multiplet, and the parameters m; = pim 
are the observed masses of the particles. The states that are described by the eigenfunctions of can be naturally 

called the mass states for any choice of the operators with continuous spectrum, that is, for any choice of the matrix 
ILtA) including the case P^ = I. 

However, we can choose a basis in the space / 2 in the form of a superposition of the mass states. Consider an 
arbitrary unitary matrix U with components U a i. The functions 


S' 


O) 


3 


(*) = UrifMjix), 


P,C,oc 


1=1 


(1.41) 


where ii x ) is the wave function of an arbitrary mass state, make up a complete orthonormal set in T-C m jy 2 




with 


the scalar product m • The basis (11.411) can be obtained as a result of the unitary transformation U of the space 


n 


(C 


4/2 


onto itself 


=u ^p,i-A x ) = I K ( x ’y)KxA y ') dy - 




(1.42) 



















The kernel of the transformation is determined by the formula 


K y) = 7A3 E E E [ d P Opiate) ® HxUy)- 

[Z7T) a=l k=l C==tl J 


(1.43) 


The elements of the basis (11.411) are no longer eigenfunctions of the operator M M determined by the mass matrix. 
The fifth operator of the complete set (denote it as J-^) is defined now as 


T (u) = UM {v) U~ 1 . 


(1.44) 


This operator can explicitly depend on the event space coordinates. It should be emphasized that the form of the 
causal Green function for Eq. (11.331) does not depend on the chosen basis 


s w(x) = _JL_ v p h f (7 ^ + W;)e i{px) d 4 P 

c [ ’ (27r) 4 fri 1 J P ' 


(1.45) 


Similarly, we can construct the representation spaces for the multiplets of the charged leptons and the down- and 
up-type quarks, where the values /q, and the matrices V, IN, M, IL for these multiplets may be different. 

Since we apply this approach to the problem of neutrino oscillations as an example, we introduce special notations 
for the multiplet of the charged leptons. The bare mass is denoted by M, without assuming in advance that M = m. 
According to the tradition in the theory of oscillations we denote the subspaces of the mass states by Greek letters. 
We will write the dilatation parameters as r]p instead of /q. Accordingly, the physical masses of the charged leptons 
will be M /3 = ppM, /3 = e, /q r. 

Applying a unitary matrix to e^ l \ we obtain new basis vectors 




= 5>< 


(e)( 1 ) 
ij C 3 ’ 


E 


= S n 


E 


a W *a i0) =6- 
u i U j ' 


j —1 *=1 f) =e,/x,r 

Therefore, the orthogonal projectors for the representation space of the charged leptons take the form 


pg> =a (/3) ®a (/3) . 


Hence, 


«“’ = E mW = E M 9 p !»- 


/9=e,/i,r ^ @=e,[i,T 

We determine the basis elements of corresponding to the basis (11.151) . as 


\P^ e ) ( x ) = %b ye l (x) a 

gX,vp y > r qX,vp y > 


(e) 


,(/?) 


where 


y qX,vp 


3/2 

f T j = _ n, ,p~ ivpfa) 


We determine the basis elements of corresponding to the basis (11.371) . as 


Kx,M 0 ( x ) = 


- „/,( e ) 


,(/3) 


(1.46) 


(1.47) 


(1.48) 


(1.49) 


(1.50) 


(1.51) 


where Mfj (x) are the plane waves describing particles with mass Mp 




(") 

pX,Mp 


(x) = 




w ( M /3) g -i (pa:). 


n pX 


(1.52) 


We do not introduce special notations for the quark multiplets, but use the appropriate indices only. 















9 


II. MODIFIED MODEL OF ELECTROWEAK INTERACTION 

Now we are able to write the Lagrangian of the modified theory describing the electroweak interaction. If, as usual, 
the interaction is carried by the gauge fields associated with the group SU(2) x 17(1), then the Lagrangian of such a 
theory is the Standard Model Lagrangian supplemented by the singlets of the right-handed neutrinos with the changes 
in the fermion sector only. 

The Lagrangian for the physical fermion fields in our model is written as follows: 

Cf = Co + Cint , ( 2 . 1 ) 


where 



E 




(0 M tf (<) )Y*!P (i) 


- ^2 

i=i/,e,u,d 


is the Lagrangian of free fields and 


( 2 . 2 ) 


Ci n t ^ (l) M (i) (i7/u)^ W 


g 

2V2 

9 

2^2 


(V (e V(l + 7 5 )+ !F (i/ V( 1 + 7 5 ) tf' (e) W+) 

0 7 m(i + 7 5) i+ ^(“) 7 ' i (l + 7 5 ) 

-eJ2 


g 


i=e,u,d 


2 COS gw . 


^ tfddy 2 ( T (i) _ 2Q« Sin 2 g w + tW 7 5 ) ^ (i) ^ 


2=is,e,u,( 


(2.3) 


is the Lagrangian of the interaction between the fermion fields and the vector boson fields W±, Z lx . A fl . and the 
Higgs field H. Here $w is the Weinberg angle, e = g sin gw is the positron electric charge, T W is the weak isospin 
(T*-") = = 1/2, T = —1/2), QW j s the electric charge of the multiplet in the units of e, and v is 

the vacuum expectation value of the Higgs field. Thus, (EB formally coincides with the Lagrangian of the Standard 
Model, supplemented with the singlets of the right-handed neutrinos (see, e.g., 0). However, the field functions 'JP') 
describe not the individual particles, but the multiplet as a whole. So it is not necessary to introduce the mixing 
matrices into £ irat explicitly. 

In the Kobayashi Maskawa formalism Q, non-diagonal matrix of the Yukawa couplings is introduced in order to 
obtain different masses of fermions and generation mixing. This procedure violates the S't/(3)-symmetry of the theory. 
The action defined by the Lagrangian of free fields (12.21) is explicitly invariant with respect to SU (3) transformations 
generated by X%,xjf\xj: u \ and X^ (see (11.241) 1. Therefore, when quantizing the model, the multiplet can be 
considered as a single particle. The one-particle states in the Fock space are defined as usual, the creation and 
annihilation operators satisfy the canonical commutation relations. However, these operators carry an additional 
discrete quantum number that is associated with the mass of the state. The multiplet can be either in one of the three 
mass states, or in a pure quantum state that is a superposition of the states with fixed masses. In a certain sense we 
may say that there exist only four fundamental fermions in our model. 

We will discuss all this in detail, using the approach described in (34|. Consider the basis (11.151) . We can write the 
components of the field functions as 


^(x) = 


Tv 


dq 3/2 (0 
- n\ 






(V s/2 tT7t7 Tv 


rfq 3/2 *( i ) 

— /V n \ 


-?,C a i,c( q ) + e lw ^) u _ gC a/- f (q) 

-9,ca7,c( < l) + e iw( ^%c*a+c(q) 


(2.4) 


With the help of Noether’s theorem we can find the integrals of motion for the fields. The analog of the energy- 
momentum tensor associated with the translational symmetry is defined by the relation 


T, 


a/3 




(2.5) 
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It ensures the existence of the conserved vector 

^/3 = E E [ 9/3^q[a+, f (q)a,; c (q)-a7 iC (q)a+(q) 
l=i c=±i J 

The current vector associated with the global gauge symmetry of the Lagrangian is defined by the relation 

3 


j a = \ x hc^i\x). 


( 2 . 6 ) 


(2.7) 


It ensures the conservation of the total field charge 

3 r r 

s = E E / dt i[ a tc(q) a i,c( t i) + 3i,c( q ) a tf( q ) 

i=i(=±r 

The tensor associated with the S'C/(3)-synnnetry of the fields is defined by the relation 
Ska = \Y,[^i V) {x)la{X^\x)). + 

A i =1 

It ensures the existence of eight integrals of motion 


k = 1...8. 


( 2 . 8 ) 


(2.9) 


X 3 = E dq a + i, C ( q ) a i,t( q ) + a i >C ( q Kc( q ) “ at,c( q ) a 2,c( q ) _ a 2, C ( q K<( q ) 

C=±i"' 


( 2 . 10 ) 


= 


C=±i ‘ 


J dq a + liC (q)a liC (q) + a i c (q)a+ c (q) + a + 2iC (q)a 2>c (q) + a 2iC (q)a+ c (q) 


- 2 H,c( q ) a 3,c( q ) - 2 a 3,c( q ) a s)t( q ) 


( 2 . 11 ) 


x ! - E 

C=±i 


^ = i E 

C=±i 

t 4 = e 

C=±i ’ 


^5 = 1 E 

C=±i 


df 6 = e 

C=±i 


*7 = i E 

C=±i 


dq 


dq 


dq 


dq 


dq 


dq 


a + u( q ) a 2)c( q ) +a~i,c( q ) a ^c( q ) + aE( q ) a u( q ) + Ef( q ) a i)c( q )_| ’ 
a + i,c( q ) a ^c( q ) + a ~i,c( q ) a ^c( q ) ~ a E( q ) a u( q ) - a_ 2,c( q ) a u( q ) > 
a + i,c( q ) a 3)c( q ) +a~i,c( q ) a Jc( q ) +aE( q ) a u( q ) + a_ 3,c( q ) a bc( q ) > 
a + i,c( q ) a ^c( q ) + a ~i,f( q ) a ^c( q ) _ a E( q ) a u( q ) _a ^,c( q ) a u( q ) > 

a gc( q )a^, c (q) +a^ 2iC (q)a+ c (q) +a+ c (q)a" x (q) +a)3 iC (q)a+ c (q) , 
[at,<( q )a^. c ( q ) + E c (q)a£ c (q) - a+ < (q)a“ c (q) - a“ 3iC (q)aJ c (q)j. 


( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 
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Using these integrals of motion and the total field charge it is possible to construct nine linear combinations 


Xik= Y 

C=±i 


dq 


a tc( q ) a fc,e( q ) + a i,c( t 0 a fc,c(*0 > 


k,l = 1,2,3, 


(2-18) 


and three of them are diagonal in the indices l, k. 

We do not consider the angular momentum tensor, because its analysis is not critical for what follows. Since the 
Lorentz group generators (ll.25[) have the standard form, then the angular momentum tensor is the same as in fdil . 

The integrals of motion for the fields are associated with the generators of the Lie algebra of the symmetry group 
of the theory as follows: 


Z&M(x) 




(2.19) 


Here Z = jp^, Q, and Z = {Vp, Q, <W....}. These relations enable one to interpret a~^(q) and a i( .(q) as 

the operators of neutrino creation and annihilation in the state l with the kinetic momentum q and the polarization 


* — 

' l,C 


(. Accordingly, a^(q) and a ( c (q) are the operators of antineutrino creation and annihilation in the state l with the 


kinetic momentum q and the polarization £. Eq. (12.1911 and the invariance condition with respect to the change of 
particles to antiparticles for Vp yield the canonical commutation relations 


h,c(q)- a fc,<'(q')J + = Si k 5 ic 6(q- q'), [a z, c (q), a£ c (q')J + = 6i k 6 a >5(q. - q'). 


( 2 . 20 ) 


Let us now consider linear combinations of these operators 

3 


a 5, c (q) = 51 t7 «; a ; ± c(q), * a l((i) = 55 


( 2 . 21 ) 


1 = 1 


1=1 


where U a i are components of an arbitrary unitary matrix LJ. The commutation relations for a^(q), a^ ^(q) are 
canonical. When expressed in terms of these operators, only Vp and Q are diagonal. However, using Q and the 
remaining integrals of motion we can always construct three linear combinations of the form 


Ka= I d <l * a+ a,c( C l) a a,c( C l)+* a 


a,£ ' 


«»C 


(q)a+ 




(q) 


a = 1, 2, 3. 


( 2 . 22 ) 


Therefore, the operators a^(q), a^ ^(q) lead to well-defined states in the Fock space. The difference between the 
number of particles and the number of antiparticles of each type a with the same kinetic momentum q is an integral 
of motion. 

Consider now the basis (11.371) . We can write the components of the field functions as 




( 2 t r ) 3 / 2 




EE 

1=1 c=±i• 

3 


dp (z) 

dp * (z) 

- n) 


( 2lr ) 3,a Si5d -M 


e Kpx) u ( p m c l) a lx {p ) + e i(px) ziL™5a+ f (p) 
3 -i (p ,)-(_^U a - c(p ) + e i&»)«^)*a+ c(p) 


The operators a,^(p), a , ^(p) arise as a result of the scaling transformation 

a z ± c (p) = (w) _ 3 / 2 a z ± c(q)’ a tf(p) = M =s/ 2 * a % c( q ) ; q = p /w- 


(2.23) 


(2.24) 


Therefore, using Eq. (12.191) and the invariance condition with respect to the change of particles to antiparticles for 
Vp, we get that these operators as well as their linear combinations 


WC 


3 

(p) = E^ a± 


1=1 


i, c 


(p) 


*± 


3 




:p) = > u< 


i=i 


*± 

J al 3 l X 


(P) 


(2.25) 


satisfy the canonical commutation relations (12.201) . 
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A similar reasoning shows that the states, which are described by the operators a^-(p), a^(p) and ^(p), a^ ^(p) 
are well-defined in the Fock space, too. However, there is an important difference. The integral of motion Vp is non¬ 
diagonal in terms of the operators a^ ? (p),a^ ^.(p). This situation is quite expected. The integral of motion Vp is not 
the canonical momentum of the field, but the “kinetic momentum”. Since the superpositions of the mass states are 
non-stationary states, the non-diagonal form of Vp reflects the fact of a possible energy transfer from one state to 
another. 

The SU (3)-symmetry of the action defined by the interaction Lagrangian (12.311 is broken even if we will admit that 
the matrices are equal. Therefore, the processes of the type ^ + 7 are always forbidden. In this case, we 

come, at least in the framework of the perturbation theory, to a model with three independent generations of fermions. 

If the matrices Y are distinct, the terms of the Lagrangian that describe the interaction via the charged currents 
immediately cause a phenomenon, which is known as the mixing of generations. The matrix of the mixing coefficients 
for quarks is an analog to the Cabibbo-Kobayashi-Maskawa (CKM) matrix 

TjCKM = v (u)t V (d) , (2.26) 

and for leptons it is an analog to the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix 

^jjPMNS = V (e)t V W. (2.27) 

The elements of these matrices can be expressed by the scalar products of the basis vectors. For example, 

^nf ) ^ a) =I/™ NS . (2.28) 

2=1 

We shall use the notation U ^) MNS = P a i for the PMNS matrix elements to omit the lengthy indices. 

In the presence of the mixing, the transition probabilities for the superpositions of the mass states cannot be 
reduced to the sum of the transition probabilities for the mass states. It should be emphasized that the total transition 
probabilities for these states can be different for the same unitary matrix U, but various complete sets of the operators 
defining the mass state. The reason is that the operators (11.441) for such states can be different. Moreover, the 
measurement results can be different for different space-time localization points of the detector. 

Just this alternative is implemented. In experiment, as it is well known, the generation mixing occurs for quarks, 
while the transfer of energy from one neutrino state to another is seen as the oscillation phenomenon. 


III. NEUTRINO OSCILLATIONS 

Now we try to apply the developed formalism to the problem of neutrino oscillations. Eq. (11.3311 is obviously 
translation invariant. Therefore, shifting the argument x^ in its solution by a constant 4-vector z A ‘, we still obtain a 
solution. For the solutions describing the mass states, this transformation leads to a trivial phase multiplier. However, 
the form of the solutions that are determined by formula (11.411) will change. That is, an argument shift generates a 
unitary transformation in the solution space. 

Let us examine this issue in detail. For clarity, we will describe the neutrino with the help of the density matrices 
of pure states. By definition, the density matrix of a pure state is g(x,y ) = <F(a;)<F(y). We assume that the neutrino 
source and the detector are at a distance L. In the area where a neutrino is produced, it is described by the density 
matrix (x,y; qX, and in the detection area it is described by the density matrix g^{x — z,y — z;q,(,&). A 
mass state can be described, for example, by a density matrix 

£> (w) (a;,y;g,C,w) = -j-L e -'(q(x-y))m ( 7 ^ + m)(l - Ct 5 7 m<S#( q)), (3.1) 

which is based on the solutions (11.151) . or by a density matrix 

g (v \x,y,P,C,mi) = ^e" lfe(l " s]) (rS l) ® n (l) ) (y M pf + mi)(l - (pi))> (3.2) 

which is based on the solutions (11.371) . In both cases there is no dependence on L. 

The situation changes in the case of superpositions of the mass states. If we take a state described by a superposition 
of the wave functions (11.151) . then the group velocities of all the mass states of the neutrino will be the same: v = q/q°. 
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So we can suppose z ^ = q^L/\q\. Therefore, the density matrix describing the neutrino state at the distance L from 
the source can be written as follows: 


£* (l,) (x, y\q, Cj Wj Mfc; a, L) 
x (p;pfc) 3/2 ( 'n (z) <8) 


1 

4g° 


" 3 

^ e -i(gx)//i+i(gy)^ fc +27riL/L (Zfc) 
_ k,l=l 


UakU 1 


( 7 /V +m)(l - C7 5 7 Vo G?))- 


Here we introduced the notation 


27r|q| 

m(mi - mfc) 


27r/3 

{mi - m fc )\/l - P 2 ’ 


P 


v . 


(3.3) 


(3.4) 


Recall that in formula (13.31) , as well as in formulas (13.11) , (13.21) , the 4-vector satisfies the condition g 2 = m 2 , 4-vector 
Sy{q) determines the direction of the polarization of the particles, £ = ±1 is the sign of the spin projection on this 
direction. If we consider the density matrix (13.31) . averaged over the parameters a and £, we obtain 


e M {x,y;q,m) 


1 

12g° 


53 e- i(9(x " y))M! /if (n w ® n (i) 
_Z = 1 


(7/V + m). 


(3.5) 


It is clear that the dependence of the density matrix on the distance from the source appears only when the source is 
coherent. This requires the dimensions of the source to be small compared with the parameters L^ lk \ 

If we take a state described by a superposition of wave functions (11.371) . then the group velocities of the mass states 
Vi = p /pi are different. However, if m 2 fc /|p| 2 -C 1, this difference is extremely small, and we can set 


z° « L(pf + p°)/(2|p|) w Lp° k / |p| « Lpt/\p\, 


So 


where 


~V -P°k) = n ik \l + 0 {ml k /\v\ 2 )), 


OjW = 


4tt|p| 


(m 2 - m 2 k ) 


(3.6) 


The density matrix describing the superposition of the neutrino mass states at the distance L from the source, can 
be calculated if one uses the explicit form of the spinors J (see EH)- With the notation (13.61) . we have 


0 (l/) {x,y,p,C> mi,m k ] oi, L) 

if 1 


-i( Pl x)+i(p k v)+2niL/LW ^(1) g, U ak U* al 


X(7 ij.Pi +mi){l - C7 5 7a 1 'S'o (pO) 


4 ^(l + 7 °) + J 4 ^(l- 7 °) 


Pi + mi 


Pi - m 


i V _L- e -i( P ^)+i(P^)+2-iL/b (fci) Ui) ®h {k) )u ak u* al 
8 V J 


Pi + mi . 0 v jPj-mi 


Pfc + m fe 


(1 + 7°) + 


P^ - m fe 


(1-7°) 


(7/iPfc + m fc)(l - C7 5 7 m*S'o (Pfc))- 


If we consider the density matrix (l3~7l) . averaged over the parameters a and £, we obtain 


(3.7) 


8 {v) V y, P, mi) = 53 TTo e 
!=1 ^ 


1 -i(p,(x-Ii)) /„(i) , 5 > * (0 


<g) (7 M pf + mi). 


(3.8) 
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Note that in the phenomenological theory of neutrino oscillations the states, which are described by the density 
matrix m are known as the approximation of equal velocities, and the states, which are described by the density 
matrix (13.71) are known as the approximation of equal momenta. Using (13.31) and (13.71) we can get the probability 
of the neutrino transition from one state to another. These states differ from each other, so for the state (13.31) the 
decoherence effect is absent and the state (13.71) spreads due to different group velocities of its components. Due to 
the fact that the interaction conserves the energy, but not the group velocity, as will be shown later, the approximate 
equality (£„ is the average value of the neutrino energy) 


l(ik) _ L (ik) _ L (ik) = 


4tt£„ 


( m i ~ ml) 

holds, and the condition of applicability of both formulas (13.31) and (13.71) are defined by the relation 


(3.9) 


m lk L 

Ip \ 2 l { osc 


(3.10) 


In accordance with the fundamental principles of quantum mechanics the probability of observing the state a at a 
distance L from the production point for the state /3 is 

= SpCeOu V, a, L)q*{x, y; 0)). (3.11) 

Therefore, in both cases 

3 

= E U pi U* pk U ak U* al e 2 ^ L ^. (3.12) 

k,l =1 


The formula itself is quite trivial. Mathematically it gives the sum of the squared absolute values of the projections 
of the new basis unit vectors in the three-dimensional vector space over the field of complex numbers onto the original 
basis vectors. In our case, the unitary transformation of the basis is generated by the phase factors arising as a result 
of space-time translations. Regardless of the initial basis, i.e. the parameters Upi , the result is the same. Therefore, 
formula (13.121) describes the evolution of an arbitrary superposition of the mass states. In particular, we can consider 
the states that are conventionally called the flavor states. Their wave functions are defined by the PMNS matrix. 

In the phenomenological theory of neutrino oscillations only the transitions between the flavor states are considered. 
That is, it is postulated that the neutrino is produced in a flavor state. This postulate is quite natural, since the PMNS 
matrix is the only non-trivial unitary 3x3 matrix that can be built using the elements of the representation space 
for leptons. However, as we have seen, the flavor states are not uniquely defined. Their properties are associated with 
the type of the mass states, from which they are composed. The only way to answer the question, in what state the 
neutrinos are produced, is to calculate the corresponding probabilities of the processes. 


IV. PION DECAY 

In the processes involving neutrinos only the charged leptons (the electron, the muon and the tau-lepton) are 
detected directly. As the charged leptons have fixed masses their description should be done using the density 

matrix of the mass states g( e \x,y;p, (,,r]p) or &( e '( x >y>P>G>Mp). These density matrices are quite similar to the 
density matrices of neutrinos that are given by m, (EH). The result of the total probability calculation does not 
depend on the type of the density matrix. It is quite obvious, since any mass state is defined by an eigenfunction of 
the operator which is determined by the mass matrix. 

In the neutral currents processes the oscillations cannot be observed, so we can sum directly over all the discrete 
quantum numbers of the neutrino final states using either the density matrix (13.51) or (13.81) (multiplied by 6) in the 
calculations. Thus, the probability does not depend on L. 

Therefore, we will examine how neutrino oscillations affect the probability of the processes which occur via the 
charged currents only. As an example, consider the pion decay 

n + =>l^ + v, l~p = fi + , e + . 

Let the 4-nromentum of the pion be k^, k 2 = ml, and 4-momentum of the lepton it be p^, p 2 = Mj. We assume 
that the distance from the source of neutrinos is L , and the linear size of the source is Lq. For clarity, we assume that 
the pion is at rest: k° = m, T , k = 0. 
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First, consider the probability of the process that produces a charged lepton with mass Mp and a neutrino with 
mass to;. The probability of this process in the Fermi approximation is given by the formula 

W « = f ^ J d<idp (4.!) 

x Sp { 0 (e ) (a;, y, -p, C, /?b M (1 + 7 5 )e M (y, x; q ,, (, l) 7^(1 + 7 5 )^^e _i(fc(a:_y)) } . 

The density matrix of the neutrino g^ v \y, x; q, £, l) can be taken either in the form (13.111 or (13.21) . The density matrix 
of the charged antilepton g^ e \x,y; — p, £,/3) can be taken either in the form analogous to (IB1I) 

Q {e) (x,y, -p,C,Vp) = ^ (a<« ® & (/3) ) ( 7 ^ - M)( 1 - C 7 5 7 'X(p)), (4.2) 

or in the form analogous to (IB2I) 

0 (e) (x, y, ~P, Ci Mp) = e <Pf>( x ~v)) fa (/3) <g>a UJ) ) ( 7 ' M p% - Mp)( 1 - C7 5 7m' S 'o (Pp))- ( 4 - 3 ) 

4 Pp v 7 

Substituting these expressions into eu and making elementary calculations, we obtain 


Wpi = ^^-P/3 iPpi Mj{rnl — Mp + mf) + mf(ml + AfJ - mf) (nil ~ M p + mf) - 4 m 2 mf. (4.4) 


This expression shows that the total probability of the pion decay, 


Wp = J2 (1 + E OiPp^ml!^)) 


(4.5) 


weakly depends on the mixing parameters Ppi (here £ v ~ ( m 2 — Mp)/2m n ). 

Now we calculate the probability of detecting a superposition of the neutrino mass states produced in the pion 
decay with the emission of a charged lepton with mass Mp at a distance L from the source. This probability is equal 
to 

W *° = j <fxd‘y J dctdp (4 6) 

x Sp {e (e) (x, 1 r, -p, C, /3)7 m (1 + 7 *)Q {v) (■ y, X] q, (, a; L) Y (1 + 7 5 )fc M fc i ,e“ } . 

First, perform the calculation taking the density matrix given by formula (13.31) as g^ v \y,x m ,q,C„a m ,L). That is, we 
assume that the group velocities of the mass states are the same. We change the integration variables 

zt=x*-y> i , zl = {x» + y»)/ 2, (4.7) 

and limit the range of integration over the variable zJ) to the size of the area in which the reaction takes place |35| . 
Then, after dividing by the volume of the reaction area, we obtain (see appendix IP 




a 2 f 2 9 

7\/f2\ 2 


87 rmi 




3 

E ■ 

k,l=l 


^iP*p k U ak U*i sin^TT Lp/Lose) 2 7Ti l/l«% 


A 3 

^Ik 


TtLq/Losc 


(l + 0(ml k /£l)). (4.8) 


Here 


A i k = (to; + TO*,)/(2 y/mim k ) > 1 (4.9) 

is the ratio of the arithmetic mean to the geometric mean of the corresponding neutrino masses. The origin of the 
factor 


sintjrLo/Lose) 

T / r d fe ) 
7TL/Q/ L/osc 


(4.10) 
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is the incomplete coherence of the source which is due to its finite size Lq <C L. After summation over the parameter 
a in formula (14.81) . there appears an expression for the total probability of the pion decay (14.51) . 

Consider the probability of detecting the neutrino flavor states at a close range L -C \Lqsc\ if the source is quite 
compact Lq <C L <C \Lqsc\- Then, replacing U a i to P a i we obtain 


W„ a = 


G\fl 


,3 “ 87TTO 3 


M p{ml ~ Mfj‘ 




Y A 3 

k,l =1 


(i + o( m f, k /£v) s ) ■ 


(4.11) 


This formula shows that not only the probability of detecting the main neutrino flavor /3 is nonzero, but the probability 
of detecting the other flavors is nonzero as well. Thus, the state of the neutrino that was produced in the decay is 
a superposition of the flavor states described by the density matrices (13.31) . It should be emphasized that we do not 
talk about a very small contribution ~ mf k /£%. 

Now take the density matrix (1371) for Q^ v \y, x\ q, £, a; L), that is, we assume that the mass states included in the 
superposition have the same canonical momenta. Using the same assumptions as in the previous case, in particular, 
the change of variables (1X71) we obtain 


wL = 


G 


2 f 2 


8-7T m 3 




t p n p; t v„ t u:, ^SS 


k,l =1 


nLo/L, 


(ik) 


(l + 0{mlJ£l)). (4.12) 


When detecting the flavor states (U a i = P a i ) at small distances from the compact source 

Wp a = ^Mj(ml - Mp) 2 (s a p + 0{ml k /£lj) . (4.13) 

Therefore, in the decay process the flavor state /3, which is described with a good accuracy by the density matrix 
(1X71) is produced. This state can be called a flavor state in the conventional sense. Naturally, its evolution is described 
by the standard formula 


V 


l'09)->''Ca) 


3 

E PpiPpkP^Ki 


k,l =1 


e 2ni L/Lg k ^ 


(4.14) 


which is analogous to formula (13.121) . 

The question, what exact states are produced in the weak decays of particles, is still open. In our model, this state 
is pure in the quantum-mechanical sense, so that when building this state we need to vary, in general, all the quantum 
numbers of the final state of the neutrino. If this state could be expressed as a finite set of plane waves, then we would 
have to choose two appropriate matrices: a Hermitian matrix IjM (see. ( 11.361 )1 and a unitarv matrix XJ^ (see. D3D). 
It is not improbable that the required state can be a wave packet, although not an abstract Gaussian wave packet, 
but a wave packet with characteristics explicitly depending on the type of the process in which it occurs. However, 
this is not essential if the condition (13.101) holds. 

When detecting the neutrino states from the very large (Lq L osc ) source we have 

/t2 3 

Wpa = MJK - M l) 2 E PfnPpiUM (l + O(L osc /L 0 ) + 0[mlJ£l)) . (4.15) 

In this case oscillations are absent. If U a \ = 5 a i we get the probability of the process that produces a charged lepton 
with mass Mp and a neutrino with mass m a . And if U a i = P a i we get the probability of the process that produces a 
charged lepton with mass Mp and a neutrino flavor state a. 

Similarly, we can consider the decay 


7T => Ip +f, Ip = n ,e , 

using the density matrices for antineutrinos (EH) , ESb EH), EH- The corresponding formulas can be obtained 

(Ik) 

by replacing the sign of the oscillation length in the expression for the probability, that is, by replacing Lose with 
Lose - As is well known, this difference in the formulas for the probabilities of the processes involving particles and 
antiparticles can indicate CP violation in the theory. 
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Thus, our approach adequately describes the phenomenon of neutrino oscillations, which is, in fact, the phenomenon 
of the oscillation of the probabilities of the processes involving neutrinos. The flavor state, constructed as a super¬ 
position of the mass states with the same momentum, is distinguished by the fact that it is very close to the state 
which is formed in the decay process. The possibility to use such flavor states in order to describe neutrinos when the 
distance from the source is much larger than the oscillation length, directly follows from the smallness of the neutrino 
masses compared to their energies. Therefore, for neutrinos, which are really observed, the phenomenological theory 
of oscillations based on the ideas of Pontecorvo (see, e.g., Q), is a very good approximation for quantum field theory. 


V. CONCLUSION 

In this paper we put forward a modification of the electroweak interaction theory, in which the fermions with the 
same electroweak quantum numbers are placed in fermion multiplets and are treated as different quantum states of 
a single particle. That is, in describing the electroweak interactions it is possible to use four fundamental fermions 
only. In this model, the mixing and oscillations of the particles arise as a direct consequence of the general principles 
of quantum field theory. This approach enables one to calculate the probabilities of the processes taking place in the 
detector at long distances from the source. Calculations of higher order processes including the computation of the 
contributions due to radiative corrections can be performed in the framework of perturbation theory using the regular 
diagram technique. 

The developed approach is used to study neutrino oscillations. It is shown by the example of the pion decay that 
the states of the neutrino produced in the decay process can be described by a superposition of states with different 
masses and identical canonical momenta with very high accuracy. 
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Appendix A: Explicit form for the generators of 5(7(3) transformations 

A direct calculation yields the explicit form of the generators of 5(7(3) transformations = KS V '> 

X[ v) = (n (1) g n (2) ) 2? ( i)2?p} + (n (2) g n (1) ) X> (2) X> ( -}, 

= i (n^ g 2?( 2 )2> ( -; - i g n (2) ) V (1 )V^y 
g g) n^ 2 ^ , 

X[ v) = g) n (3) ) 2? ( i)2?p} + (n^ g n (1) ) V (3) V~^ 

X£ ] = i (n (3) g n (1) ) 2?(3)X> ( -5 - i (?r (1) g n (3) ) V {1) D^ y 
X^ = (n^ g ^ V (2) V( 3 ] + (n^ g n (2) ) V {3) V~], 

X^ = i (n (3) g n (2) ) V (3) V^ - i (ra (2) g n (3) ) V (2) T>q], 

Ag") = g g n^ 2 ^ — 2 g . 
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Operators of finite transformations are given by 


where 


d'V) = 

^n® 0 n® 

+ 

^n® 0 n® 

+ 

n ® 0 n®^ 

II 

^n® 0 n® 

+ 

^n® 0 n® 

+ 

n ® 0 n®^ 


• v(y ) • 

cos a — iX i sin a, 

• vi 1 ') • 

cos a — 1 X 2 sin a, 


u. 


O) 


(a) = 0 h 

U^\a) = (n^ 0 n 


(a) = ( 

3 W («) = ( r 


= In® 0 n 


(3) 

*( 2 ) 

( 2 ) 


^ + (n® 0 h' ' | e 


Uk”’(a) = I n® 0 n® 


U^”\a) = ^n® 0 n® 


*U)j „—ic 
*(!)' 


^n® 0 h 


n®) e ic 


*(3) 


^n® 0 n® j + ^n® 0 n 
^n® 0 n®^ + ^n® 0 n®^ 

"® j ( n ^ ® 

0 n®^ 


^n® 0 h 
^n® 0 n® ^ 


• v'( I/ ) • 
cos a — 1 A 4 sin a, 

• v( u ) • 

cos a — 1 A 5 ' sin a, 

• \r(is) • 

cos a — iAg sin a, 

• vM • 

cos a — 1 A 7 sin a, 


U&\a) = (n® 


0 n® ) e 2ia 


^n® 0 n®^ + | 'n ® 0 n®^ 


(n (1) 0n (1) ) = ^(21 + X® + 3X®), (n (2) 0B (2) ) = i(21 + X® - 3X ®), 

^n® 0 n® j i (l — A ®). 


(A2) 


(A3) 


Appendix B: Density matrix of antineutrino 


(Bl) 


The mass states of antineutrino are described either by the density matrix 

f?® (x, y\ —q, Cj /m) = (n (l) 0n (l) j ( 7 ^ - rn)(l - C7 5 7a»-So (g)), 

based on the solutions similar to (11.151) . or the density matrix 

£® (x, y, —p, C> mi) = e '( p i( x ~ y )) ^ n (0 (g, j (7 M pf-mj)(l-C7 5 7 ^(pj))- ( B2 ) 


based on the solutions similar to (11.371) . In both cases, there is no dependence on L. 

If we consider a state, which is described by a superposition of wave functions similar to (11.151) . then the density 
matrix for the state of an antineutrino at a distance L from the source can be written as 


1 


e w {x,y,-q,£,nuHk',(x,L) = 4 ^ 0 


x (niHk) 3/2 (n (l) 0 n®) U ak U* 


CJ i(qx)m-i(qy)/j. k +2mL/L , -‘ k '> 


E 

k.l=l 


{ln<f - rn)( 1 - C7 5 7 /x<S'o (?))■ 


(B3) 


If we consider a state, which is described by a superposition of wave functions similar to (11.371) . then the density 
matrix for the state of an antineutrino at a distance L from the source can be written as 

0 ® (x,y;-p,(, mi ,m k ;a,L ) 


1 


= - y 


1 


0 i(pix)-i(pi c y)+2mL/L (lk '> 


8 1 

x(7^r - m 0(! - C7 5 7 p'S'o (pO) 

3 


( 8 ) 




«z~(i +7 » )+ /4±~r, ..0 

Pi - mi K \ Pi + m i 


u ak u* al 

(l — 7°) 


(B4) 


e 1e 


1 


3 -i(pix)+i(p fc i/)+ 27 riL/L ( ' ci) 


8 *1=1 vtM 


0 n®^ 








(7MPfc - m fe)( 1 - C7 5 7 /x*S'o (Pfc))- 
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Appendix C: Calculation of pion decay probability 


The expression (14.61) after summation over the polarizations of the final particles can be easily reduced to the 
following form 

wL = 


% = 4 ^ J J r fc o J d 4 xd 4 y J d 4 qd 4 pS(p 2 - M%)6(q 2 - to 2 ) 


e i(qx)[m k /rn]-i(qy)[rni/rn]+i((p-k)(x-y))+2'rriL/L ( ' lk) 

k,l =1 

x ( mi/m) 3/2 ( m k /mf !2 PpiPp k U a kU* t 

x Sp [( 7 ^ - M /3 ) 7p ( 1 + 7 5 )( 7 , 9 " + m) 7 A(l + 7 5 )k p k x ] . (Cl) 
Having calculated the trace of 7 -matrices and changed the integration variables 

q ‘ 4 -A q^(m/rriki), m k i = (m k + mi)/ 2, (C2) 

we get 


/3 “ (2tt) 6 /c° 


E 


p P ip; k u ak u* al 


A 3 

UM=i 


J d 4 qd 4 p5(p 2 - Mp)5(q 2 - m 2 kl ) 


X (2 ( P k)(qk) - k 2 { P q)) J rf4 a . d 4 2/e i(^)[m fc /m fci ]-ito)[m 1 /™u]+i((p-fe)(^-y))+2 7 rii/i (ifc) 


Here 




47r| q| TOfc+TO; 

A/fe = 




We change the integration variables 


(*°,x) 
.0 


2 ^/mim k 

a /0 / /\ / /D U(UX ) 

u°x ° + (ux ), X + UI ° + 


■ (C3) 


(C4) 


1 + u° J 


where 


(i/°, y) ( A'° + (uy'), y' + uy' 0 + ) > 

V 1 + rt u / 


w^ = q^/rriki- 


(C5) 


(C 6 ) 


The Jacobian of transformation (IC5I) J = 1. This change of variables yields the following expression in the exponent 

(( * - pM - m 4 ) - (iitzM -m,V-(Yk-p- ~ p)g) + - p °) q ) (X' - y')) . (C7) 


/ V TO fcZ 

We change the variables 


m k i {m k i + 9 0 ) 


z'_=x' —y' ; z' + = (x' + y')/2, z'° = x'°-y' 0 , z% = (x 10 + y'°)/2. 


(C 8 ) 


z! = ( I - „, q ® q ) z". 


g°(m fe j + <?°) 

The Jacobian of transformation (IC9I) J' = m k i/q°- As a result, the expression in the exponent takes the form 


(( k-p)q) 


m. k i 


- m k i z'° - (m k - mi) z+ - ( k - p - 


P \ // 
z 


(C9) 


(CIO) 
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To get the probability of the process in a unit volume per unit time we use the standard Fermi ansatz. We integrate 
(1C3I) over d 4 z' + in a finite four-dimensional region V'T' and divide by the invariant volume V'T' = VT. As a result 
of this operation, a factor 


R = 


sin ((mi - m k )T'/ 2) 
(mi-mk)T '/2 


arises. It follows from Eq. (IC'51) that 


Assuming L = vT, v = q/g°, we get 


T ,_ q° T - (qL) 

mid 


rpi m ki rp _ rriki T 

1 ~ — (T 1 = TT l “’ 

q° |q| 


where Lq is the linear dimension of the region where the pion decays. Hence, 

_ sm(TrL 0 /LW) 
ttLq/L^) 

The integration over dz'_ yields the (5-functions associated with the momentum conservation 


(27t) 3 (5 ( k — p 


k°~p° 


It follows from Eq. (1C15II that 


and we have 


-q • 


, k° - p° 

P = k -^o—q, 

_ mu = *i (t » 
m k i q v 


Therefore, using an integration variable 


z°_ = z'°^ = z'°J' 
<7 U 


and integrating over it, we obtain the 5-function associated with energy conservation 

27 t5 (fc° — p° — q °) . 

Thus, 

3 PpiP^UM 

.k,l =1 


(Cll) 

(C12) 

(C13) 

(C14) 

(CIS) 

(C16) 

(C17) 

(CIS) 

(C19) 


f>a (2tt) 2 /c 0 


E -*- P iJ - Bk^ a i I 4 4 

- a "3 - / d qdp 


A 3 

^Ik 


sin(7r Lq/LW) 2niL/L m 


ttLo/LW 

x S A {k — p — q) S(p 2 - Mj) 5{q 2 - m 2 kl ) (2 {pk){qk) - k 2 (pq)) 


(C20) 


The integral over the momentum variables coincides with the one that appears when calculating the probability of 
the mass state creation. Therefore, 


r 2 f 2 

WL = fU 


0a (2t r) 2 k° 


Ppi P pk U »kU*i r 4 sin(71 -L 0 /LW) 27liL/L m 2 2 \ 

^ Af k J dq nL 0 /Lm 6{q mu) 

x 6(m 2 — Mj + m 2 kl — 2 (kq)) [Mp(m 2 — Mj + mh) + mh(m 2 + — m 2 kl )\ 


(C21) 


Assuming k = 0 and neglecting m k i in comparison with the neutrino energy we obtain the formula (14.81) . 
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It should be noted that, the change of variables (TOT) and limiting the integration domain to the size of the area 
where the reaction occurs is not a clearly defined operation, but rather a prescription. However, such_a procedure is 
not only justified by physical considerations of Schwinger f35j, but is also mathematically consistent [HI]. 
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